Characterizing the structure of small-world networks 
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We give exact relations which are valid for small-world networks (SWN's) with a general 'degree 
distribution', i.e the distribution of nearest-neighbor connections. For the original SWN model, 
we illustrate how these exact relations can be used to obtain approximations for the corresponding 
basic probability distribution. In the limit of large system sizes and small disorder, we use numerical 
studies to obtain a functional fit for this distribution. Finally, we obtain the scaling properties for 
the mean-square displacement of a random walker, which are determined by the scaling behavior of 
the underlying SWN. 

PACS numbers: 05.10.-a, 05.40.-a, 05.50.+q, 87.18.Sn 
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Networks occur in many contexts in the sciences and 
humanities. Among these are neural networks social 
networks food webs and computer networks 

|^-^| . Many properties of such complex systems can be 
understood by considering the network of interactions 
which connects their components. For this reason, much 
recent research has focused on the structure of these net- 
works . 

The study of network models has recently been ad- 
vanced by the introduction of small-world networks 
(SWN's) by Watts and Strogatz @]. These networks 
have an ordered structure locally, but are random on a 
global scale. This combination of features suggests that 
SWN's can be used to describe the behavior of many 
real- world interacting networks [ fL5| . However, as origi- 
nally defined, SWN's do not exhibit a property found in 
some large-scale networks such as the world-wide web, 
in which the probability distribution for the number of 
nearest-neighbor connections to a given site (the 'degree 
distribution') is a power law, rather than Poisson-like, as 
in the SWN's (l(J. Alternative models which do exhibit 
this power-law distribution have been proposed [|lG] . 

In this article, we study SWN's where shortcuts ac- 
cording to a given degree distribution are superimposed 
on a regular network. For this class of networks, we show 
some exact results which are independent of the specific 
degree distribution employed. These can be used to find 
approximations for the basic probability distribution for 
the network. Furthermore, on the basis of numerical sim- 
ulations of the original SWN model (with Poisson-like 
degree distribution) , we propose that the basic probabil- 
ity distribution is characterized by a function of a sin- 
gle variable in the limit of small disorder (as defined be- 
low) and large system size: the so-called 'small-world' 
regime. Using this functional form, various structural 
properties of networks in the small-world regime can be 
accurately determined without the need for further, ex- 
pensive computer simulations. We also study how the 
structural properties of this particular SWN model af- 
fects an important example of dynamics on SWN's: the 



time-dependent mean-square displacement of a random 
walker. 

The SWN's are defined by starting from a one- 
dimensional regular network with periodic boundary con- 
ditions and L — 2N nodes, each node being connected to 
its 2k nearest neighbors. We add shortcut ends to each 
node according to a given degree distribution by follow- 
ing the prescription of Ref. jTsj. For the original SWN 
model, the degree distribution is V q — (1 — p) V q (kp) + 
p 1 D ( q -i){kp), where V q {\) — exp(— \)\ q /q\ is the Pois- 
son distribution. We then select pairs of shortcut ends 
at random and connect them to each other, thus creat- 
ing a shortcut. This procedure (with the above T> q ) of 
generating the networks is equivalent to the procedure 
outlined by Newman and Watts M , with kp as the prob- 
ability of having a shortcut at a given site. Thus, on 
average, there will be x = kpL shortcuts in the network. 
The present numerical results were calculated for k = 1; 
however, generalization to other k is straightforward and 
indicated in our equations (which we have verified nu- 
merically). The 'small disorder' regime corresponds to 
p <C 1, and we will work in the limit of large system sizes 
such that terms of 0(1/ L) can be omitted. 

The SWN is characterized by two types of distances: 
the 'Euclidean distance,' defined as the shortest distance 
between two sites before any shortcuts are introduced 
into the network; and the 'minimal distance,' which de- 
notes the shortest distance after shortcuts have been 
added. An important probability distribution in SWN's 
is then P(m\n), defined as the probability that two sites 
have minimal separation m, given that their Euclidean 
distance is n. As shown previously C3), P(m\n) can be 
written in the form 



P(m\n) 



f(m), 



m < n, 



1 - EIU 1 /(*), m 



0, 



n, 
m > n. 



(1) 



Thus, P(m\n) is fully characterized by the function f(m) 
p0| which can be regarded as the probability that two 
diametrically opposite sites have minimal separation m. 
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The above result is exact for SWN's with arbitrary degree 
distributions. 

Many structural properties, such as the average min- 
imal separation between two randomly chosen points, £, 
can be expressed in terms of /(to) H], and these expres- 
sions can be used to derive exact relations connecting var- 
ious quantities of interest. An important relation which 
can be derived in this manner, relates the average num- 
ber of sites within m hops of a randomly selected origin, 
V(m), to F(m) = J2T=i /W> * ne cumulative probability 
distribution, as 



V (to) = 2km +[L- 2km]F(m), 



(2) 



Through our simulations, we have confirmed the validity 
of Eqs. (0) and @ for several different types of degree 
distributions, including networks with power-law degree 
distributions. These exact relations can be used to ob- 
tain quantities of interest for the network. In a simple 
model for disease-spreading |l2] |, for example, V(m) cor- 
responds to the number of infected sites after m time 
steps. Furthermore, Eq. (|J) tells us how to calculate 
f (to) [and hence /(to)], once we know V(m). Thus, any 
approximation for V(m) can be converted into an ap- 
proximation for /(to). In the remainder of this article, 
we will illustrate this procedure by focusing on the orig- 
inal SWN model. 

We will now use the mean-field result of Newman et. 
al for V(to) to obtain an approximation for /(to), and 
hence, to motivate a new scaling form for this function. 



In the asymptotic limit (x 3> 1), using Eq 
tion with the mean-field results of Ref. [13 
relation between /(to) and F(m): 



1) in conjunc- 
, we obtain a 



/mft(TO) = 4 k 2 p F mft (m) (1 - F m{t (m)), 



(3) 



which is the logistic growth equation, as pointed out 
by Barbour et al pi] . Eq. (^) can easily be solved 
for f m {t(m) f22| . We note that the final expression for 
/mft(TO) is consistent with the scaling form /(to; N, p) = 
h (to /N, pN) /N derived in Ref. @ . 

We now consider the central moments of this distribu- 
tion, defined as 



N 



(to)] /mft(m). 



(4) 



The scaling form for /(to) implies that the scaled mo- 
ments ni/N are functions of x only; however, f m ft(m) 
has an additional property: its central moments in this 
case take the form 
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where eg is the £ th central moment of the logistic distribu- 
tion. We have confirmed through numerical simulations 
that Eq. (||) holds for at least the first four central mo- 
ments of /(to), however, with coefficients eg which differ 
significantly from those obtained from the logistic growth 
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FIG. 1. Plot of the universal function g(y) [cf. Eq. (^J)] for 
three choices of system size L and probability of shortcuts, 
p (see text): (i) L = 10 7 , p = 10" 4 (full line); (ii) L = 10 7 , 
p = 10" 3 (stars and full line); and (iii) L = 2 x 10 6 , p = 1CT 3 
(circles and full line). 

equation. For example, from our simulations, we find 
C2 ~ 1.21, whereas the corresponding mean-field value 
is c™ ft w 0.82. Eq. (H) also implies that fig ~ p~ l , and 
hence, for fixed p, is independent of N. Consequently, 
the probability distributions /(to), for fixed p and dif- 
ferent values of N, are simply related by a translation 
of the mean values. We have confirmed this translation 
property by simulations carried out for a wide range of 
parameters in the small-world regime. 

On the basis of these features, we propose the following 
scaling form for /(to) in the asymptotic limit (i> 1): 



f(m;N,p) = g(^r(m- (m)) 



(0) 



where g(y) is a universal function. It is readily seen that 
the above scaling form is consistent with Eqn. In 
Fig. [I], we show calculations of /(to) for three differ- 
ent combinations of the parameters (p,L): (10 _4 ,10 7 ), 
(1(T 3 , 10 7 ), and (1CT 3 , 10 6 ). For all three choices of pa- 
rameters, /(to; N,p) can be collapsed onto a single curve. 
Thus, numerical simulations for just one choice of the 
parameters (p,L) suffice to determine g(y). To fully de- 
termine /(to) we need to find functional forms for both 
g(y) and (to), which we now proceed to do. 

We can represent g(y) using a two-parameter mixture 
of logistic distributions: 



(5) g(y) 



2(1 + 7) 



sech 2 (y + 7/?) + 7 sech 2 (y — 0) 



(7) 



This form ensures that g{y) is normalized, has zero mean, 
and has the correct behavior for large \y\. In Fig. |^, we 
show a comparison of the mean-field result, the fitting 
function, [Eq. (§) with 7 = 0.492 and (3 = 0.782], and 
the numerically obtained g(y). Evidently, this function 
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FIG. 2. Comparison of the universal function g(y), as ob- 
tained from numerical calculations (full line), mean-field the- 
ory (dashed line) , and our two-parameter fitting function Eq. 
(M) (stars and full line). 

describes the numerical g(y) very well, while the mean- 
field g(y) agrees less well. 

To determine (to) , we use the exact relation, which was 
proved in |l4| ] and is valid for a general degree distribu- 
tion, 



1 



(8) 



together with the definition /i2 = (to 2 ) — (rn) 2 . Thus, 
once we have fitting forms for I and /12 , we can determine 
(to). From our simulations, we have already determined 
/j,2 (see above). As for £, we find numerically that 



k£ _ log(2x) 
T ~ Ax 



a 

■ 

x 



(9) 



with a = 0.13. This numerical estimate agrees well with 
the mathematically rigorous result derived in Ref. pj| . 
In Fig. ||, we show a comparison between the fitting form 
@ , the numerically obtained values for £, and the mean 
field result, k£ m f t /L — > log (2^) /(4a;). As the inset shows, 
Eq. (||) gives £ with a relative error of less than 1%. 
Keeping only terms through 0(l/x), and recalling that 
H2IL 2 ~ 1/a; 2 , we find that (to) « I. Hence, (to) is given 
by 



(to) 



L 
~k 



\og{2x) + 0.13 

Ax x 



(10) 



The three equations (§), (0), and (|l(]) provide a full 
functional form for /(m) in the asymptotic scaling regime 
for a SWN with degree distribution T> q . Note the striking 
pattern of successive simplifications during our analysis 
of the probability functions. Initially, P(m\n) depended 
on the four variables m, n, p, and N. The exact results 
obtained in Ref. MM reduced this dependence to a 



CL 



V 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 











! % 








- 1 




; * 





10 10 10 




* numerical data 

- ln(2x)/4 

- ln(2x)/4 + 0.13 



10" 



10 1 



10 J 



10" 



10' 
x=p*L 

FIG. 3. Numerically computed average minimal distance 
between two randomly chosen nodes, i, (stars) using p- values 
ranging from 10 -4 to 10 -3 and system sizes L from 10 4 to 
10 7 . Dashed line: mean-field result £ m n/L = log(2a;)/(4a;). 
Solid line: our fitting form £/L — imtt/L + 0.13/a;. The inset 
shows the absolute relative error between the numerical data 
and our fitting form. 

function of three variables /(to; N, p). Furthermore, the 
scaling properties discussed in |14[| imply that, for p -c 1, 
/(to; N,p) = h(m/N,pN)/N, a function of only two vari- 
ables. Finally, we have shown in the present work that 
in the limit x 3> 1, the small- world network is charac- 
terized by a function of a single variable, g(y). This 
functional form can be used to obtain detailed, accurate, 
structural information about networks in the small-world 
regime without having to perform computationally ex- 
pensive simulations. 

The scaling forms for f(m;N,p) can also be used to 
gain insight into dynamics on SWN's. As an example, 
we consider the dynamics of a random walker on SWN's, 
a problem which has been studied by several authors 
|p3| psf . Here, we focus on the mean-square displace- 
ment, Jr 2 } of a random walker on a SWN, as a function of 
time |2^|. At long times, (r 2 ) must saturate to the finite 
value £ 2 (L,p), since, in the long-time limit, each node of 
the network has equal probability of being occupied by 
the random walker. The other length scale determining 
the behavior of the random walker is £ = 1/p, the aver- 
age distance the walker travels to reach a shortcut. Note 
that in the scaling limit, £ 2 (L,p)/L 2 , and consequently 
S, 2 / £ 2 {L,p), are functions of x only. Correspondingly we 
can write down the scaling ansatz: 



fill 



Because of the diffusive behavior for small y and satura- 
tion for large y, we expect R(y) = y for y <C £, 2 /£ 2 , and 
R(y) = 1 for y ^> 1. We have numerically confirmed the 
scaling collapse implied by Eq. (pd|) for a wide range of 
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FIG. 4. The mean-square displacement, (r 2 ), for a ran- 
dom walker on SWN's with x = 100 (upper curve), using 
L = 5 x 10 3 (circles and full line) , and 1 x 10 4 (stars and full 
line). By contrast we also show (r 2 ) for a regular network 
(x = 0; dashed line). 

x values, and an example is shown in Fig. ||. Note that we 
can readily calculate the length scale £ 2 using our func- 
tional form for f(m), illustrating how knowledge of the 
underlying structure helps in understanding the dynam- 
ical properties. 

In summary, we have given exact results which are 
valid for SWN's with an arbitrary degree distribution 
superimposed on a regular lattice. We have also shown 
that in the small- world limit (x ^> 1), the structure of the 
original SWN's can be characterized by a single function, 
g(y), and the mean value, (m), and we provide empirical 
fitting forms for both quantities. In addition, we have 
shown that the mean-square displacement of a random 
walker exhibits scaling properties which follow from those 
of the underlying SWN. The structural aspects of SWN's 
discussed in this work should be useful in understanding 
many other properties of dynamics on SWN's. 
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